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Abstract 

We present explicit expressions of the helicity conservation in ne- 
matic liquid crystal flows, for both the Ericksen-Leslie and Landau-de 
Gennes theories. This is done by using a minimal coupling argument 
that leads to an Euler-like equation for a modified vorticity involving 
both velocity and structure fields (e.g. director and alignment ten- 
sor). This equation for the modified vorticity shares many relevant 
properties with ideal fluid dynamics and it allows for vortex filament 
configurations as well as point vortices in 2D. We extend all these re- 
sults to particles of arbitrary shape by considering systems with fully 
broken rotational symmetry. 



1 Introduction 

Several studies on nematic liquid crystal flows have shown high velocity gra- 
dients and led to the conclusion that the coupling between the velocity it(x, t) 
and structure fields is a fundamental feature of liquid crystal dynamics [27]. 
This conclusion has been reached from different viewpoints and by using dif- 
ferent theories, such as the celebrated Ericksen-Leslie (EL) and the Landau- 
de Gennes (LdG) theories [4, 32, 31, 2]. Evidence of high velocity gradients 
also emerged [24] by using kinetic approaches based on the Doi model [10]. 
The essential difference between EL and LdG theories resides in the choice of 
the order parameter: while EL theory for rod-like molecules considers the dy- 
namics of the director field n(x, t) and it is successful in the description of low 
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molar-mass nematics, the LdG theory generahzes to variable molecule shapes 
by considering a traceless symmetric tensor field Q(x, t). In the presence of 
high disclination densities the LdG theory is more reliable, since molecules 
may easily undergo phase transitions (e.g. from uniaxial to biaxial order) 
that are naturally incorporated in the theory. However, the dynamical LdG 
theory is not completely established and different versions are available in 
the literature [1, 28, 29]. Here, we shall adopt the formulation of Qian and 
Sheng [29], which will be simply referred to as LdG theory. 

This paper considers both EL and LdG theories and it shows how the 
strong interplay between velocity and order parameter field refiects naturally 
in the helicity conservation for nematics. In this paper, the term "helicity" 
stands for the hydrodynamic helicity and not the helicity of the single liquid 
crystal molecule. The helicity conservation for incompressible liquid crystal 
fiows arises from the simple velocity transformation u C = u + J , where 
the vector J depends only on the order parameter field. The covariant vector 
C is the total circulation (momentum per unit mass of fiuid) and J is the 
circulation associated with entrainment of fiuid due to its local interaction 
with the nematic order parameter field. We shall show how this change of 
velocity variable takes the equation for the ordinary vorticity a; = V x w 
into an Euler-like equation for the modified vorticity oJ = V x C, thereby 
extending many properties of ordinary ideal fluids to nematic liquid crystals. 
The helicity is then given by 

n = J {-lj-c) d^x, 

and this quantity naturally extends the usual expression J (u) ■ u) d^x for the 
helicity of three dimensional ideal fiows. We recall that in Hamiltonian fiuid 
dynamics the conservation of the hydrodynamic helicity is strictly associated 
to the Hamiltonian structure of the equations and holds for any Hamiltonian. 
This point will be further developed in the last part of the paper, where all 
the results will be derived directly from the Hamiltonian structure of the 
liquid crystal equations, see [7]. Invariant functions like the helicity are 
called Casimir and are of fundamental importance for the study of nonlinear 
stability. For two dimensional fiows, such invariant quantities are given by 

I HuJ) d\ (1) 

where $ is an arbitrary smooth function. As we shall see, the same circulation 
concept leading to hydrodynamic helicity applies quite generally in complex 
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fluid theory and is related to an analogy between complex fluids and non- 
Abelian Yang-Mills fluid plasmas [8, 17]. 

In addition to helicity conservation, we present the existence of vortex- 
like configurations for the modified vorticity TD. Vortex structures are well 
known to arise in superfluid flows and their behavior is often reminiscent of 
disclination lines in liquid crystals. However, here we shall consider vortices 
that are characterized by a combination of velocity and structures fields. 
After extending these results to fluids with molecules of arbitrary shapes, 
the end of this paper discusses the geometric basis of the present treatment. 



2 Director formulation 

In the context of EL theory, disclinations are singularities of the director field 
and thus their dynamics is related to the evolution of the gradient Vn. This 
relation has been encoded by Eringen [11] in the wryness tensor 

7el = n X Vn , or (7^0" = e'^'^'rib x diU, , (2) 

which identifies the amount by which the director field rotates under an 
infinitesimal displacement dx. Thus, the EL wryness tensor 7el determines 
the spatial rotational strain [14]. In this paper we shall investigate the role 
of the EL wryness tensor in helicity conservation and vorticity dynamics 
in the EL theory. For this purpose, we ignore dissipation and concentrate 
on nonlinearity. This simplifies the resulting formulas. We also restrict to 
incompressible flows to ignore ordinary fluid thermodynamics. 

Upon denoting by J the microinertia constant [11], we introduce the 
angular momentum variable 

ctel = Jn X Dtn (3) 

that is associated to the director precession. While Dtn = dtn+u Vn denotes 
material time-derivative, the spatial derivatives of the director field will be 
denoted equivalently by dn/dx^ = diU = ri j depending on convenience. 
Upon using Einstein's summation convention, one can express the Ericksen- 
Leslie equations as [14, 6, 7] 

(dF \ 
Vn^-- — -Vj9, V-u = (4) 

c^tCTEL + -Vctel = h X n, 9tn + w ■ Vn = J^Vel X n (5) 
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where p is the hydro dynamic pressure by which incompressibihty is imposed 
and h is the thermodynamics derivative representing the First Law response 
in energy to changes in the director field 

dF ^ f dF 
h := a 



dn V (9n 

The quantity F is taken to be the Oseen-Zocher-Frank free energy 

F = Ki(divn)2 + K2{n ■ V x n)^ + K^ln x V x . (6) 

Of course, this choice is not a hmitation, because our considerations apply 
to a generic form of F. For example, effects of external electric and magnetic 
fields may be taken into account with easy modifications. The Ericksen-Leslie 
fluid equations follow immediately from equations (4) and (5), as shown in 
[6]. 

In Ericksen-Leslie nematodynamics, the quantity ctel-Tel denotes the vec- 
tor of momentum per unit mass, with components (cTel ■ 7el) j = ('^el)^ (7el)"- 
We consider the vector Cel defined as the sum 

Cel ■■=u + a^i^ ■ 7el = w + ctel ■ n X Vn , (7) 

reminiscent of the minimal coupling formula in electromagnetic gauge theory. 
We observe the following equation of motion [7] (see appendix A.l): ^ 

dtC^^-uxV xC^^ = -V{^ + u-C^^) , (8) 

where 

<l> = P + F-^\u\'-^\cT^,\' . (9) 



^ The same idea has been apphed in superfluid plasmas, that is, in superfluid so- 
lutions whose charged condensates are coupled electromagnetically [19]. However, the 
similarity with gauge theory does not end with the electromagnetic analogy. The equation 
corresponding to (8) also follows by inspection for a Yang-Mills fluid plasma (chromohy- 
drodynamics, cf. [8]) either from equation (2.35) or (2.49) of [17]. By this observation, 
chromoliydrodynamics acquires a circulation theorem and the theory of complex fluids 
inherits an analogy with Yang-Mills fluid plasma, first noticed in [14]. This minimal cou- 
pling argument requires the wryness tensor 7el to be a connection one-form: although this 
is not the case for the expression n x Vn, a connection one-form can be obtained by the 
addition of terms parallel to n. By good fortune, these extra terms make no contribution 
in (7) because ctel ■ n = 0. 
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At this point, taking the curl of equation (8) yields the Euler-like equation 



StL^EL + V X (u X CJel) = (10) 

for the modified vorticity 

Uel := V X Cbl = a; + V X (ctel ■ 7el) ■ (H) 
Notice that the velocity u can be expressed as 

U = - V X V = -V X A^^oJel + CTel ■ 7el + V(/9 , (12) 

where V = A~^ci; denotes the velocity potential, which is given by the con- 
volution of the vorticity uj with the Green's function of the Laplace operator 
(analogously for A'^uJel)- Here the pressure-like quantity is a scalar func- 
tion arising from the term V x V x A^^ (ctel ■ 7el) = ctel ■ 7el + V(/? and 
whose only role is to keep the velocity u divergence free, so that V ■ (ctel ■ 7el) 
—A(p. The relation (12) can be inserted into equations (5) so to express 
the EL equations in terms of the modified vorticity oJel- An explicit ex- 
pression of the quantity ctel ■ 7el arises from the definitions (2) and (3): 
CTel ■ 7el = </ Vn ■ An. 

At this point we recognize that equation (10) possesses all the usual prop- 
erties of Euler's equation. For example, Ertel's commuting relation 

[A , ZI^EL • V] a = A (^ ■ V«) - ZJ ■ V( A «) = (13) 

follows easily by direct verification, for any scalar function a(x, t). Moreover, 
one has the following Kelvin circulation theorem [7] 



^ (h (it + CTel ■ 7el) " dx = , (14) 

at ./r(t) 



where the line integral is calculated on a loop r[t) moving with velocity u. 
Also, conservation of the helicity [7] 



"Hel = J <^EL ■ Cel d X = y (u + CTel ■ 7el) ■ V X (w + CTel ■ 7el) d X 

follows from the relation 

dt (Cel ■ ^^el) + V ■ ( (Cel ■ ^el) u) = - cJel ■ V0 , (15) 
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which is obtained by using equations (10) and (8). Integrating equation (15) 
over the fluid volume yields 

^^EL = - ^ cJel ■ dS - ^ (Cel ■ cUel) u ■ dS (16) 

where S is the surface determined by the fluid boundary. Consequently, the 
right hand side of equation (16) vanishes when oJel and u are both tangent 
to the boundary, thereby producing conservation of "Hel- Remarkably, one 
can show that helicity conservation persists for any free energy F(n, Vn), 
that is the helicity T^el is a Casimir for EL dynamics, see [7]. At this point, 
a question about boundary conditions arises: while the condition of veloc- 
ity tangent to the boundary is the usual condition in hydrodynamics, the 
condition 

V X (n + ctel ■ n X Vn) ■ dS = (17) 

emerges here for the flrst time. Upon denoting tt = ctel x n, one has 
V X (it + ctel ■ n X Vn) = to + VtTq x V?2a, so that the boundary condition 
(17) reads as Vvr^ x Vria ■ dS = — a; ■ dS. This relation evidently differs from 
the usual "anchoring" boundary conditions (see e.g. [30]) that are widely 
used in the literature and for which the director alignment at the surface is 
insensitive to the flow. Indeed, the physical relevance of the boundary con- 
dition (17) resides in the fact that it involves both fluid and fleld variables, 
contrarily to other commonly available boundary conditions. The complete 
physical justiflcation of (17), however, requires more study in the future. 

One of the most relevant consequences of equation (10) is the existence 
of singular vortex-like conflgurations in Ericksen-Leslie nematodynamics. In 
two dimensions, equation (10) has the usual point vortex solution 

N 

ujEL{x,y,t) = ^Wi5{x - Xi{t)) 5{y - Yi{t)) , 
1=1 

where [Xi, Yj) are canonically conjugate variables with respect to the Hamil- 
tonian = ^(A~^a;)(Xj, Fj). Here ip is the potential of the velocity u, 
satisfying EL equation (4). Upon using relation (11) and suppressing the EL 
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label for convenience, one expresses the Hamiltonian as 
^(X,, Y^,CT,n)=-^J2ij2'^^^ ^'^ log I (^/^ - ^^^^ - Yk) I 

h \ k 

+ Whj {a,, 7'^} (x, y') log | (x' - X;,, y' - Y,) \ dx' dy'^ , 

where \ {x, y) \ = a/x^ + and {■, ■} denotes the canonical Poisson bracket in 
(x, y) coordinates, arising from the 2D relation V x (cr ■ 7) = V x ( Vn ■ tt) = 
{tTq, Ua], where tt = cr x n. 

Other vortex- like structures are also allowed by equation (10), e.g. vortex 
cores and patches. In three dimensions, the vortex filament 

u;,,(x,t) = I 5(x-R(s,t)) ds 

is also a solution of (10), with dtH = u(R, t). The existence of these vortex 
structures (including vortex sheets) rise natural stability questions concern- 
ing possible equilibrium vortex configurations. Instead of pursuing this di- 
rection, which will be the subject of our future work, the next sections will 
show how all the above observations also hold in the LdG theory and for 
fluid molecules of arbitrary shape. 

To conclude this section, we emphasize that in all the above discussion 
the velocity and the structure fields are strongly coupled together. Indeed, 
the singular vortex structures only exist for the vorticity cj^l — V x 
while there is no way for the ordinary vorticity uj = V x u or the 'director 
vorticity' V x (ctel ■ 7el) to be singular. This strong interplay between the 
macro- and micromotion is the same that emerges in many of the experiments 
and simulations reviewed in [27]. 



3 The alignment tensor 

In the preceding section, we investigated the hydrodynamics of a uniaxial 
nematic liquid crystal. At this point, it is natural to argue that in the 
presence of disclinations the molecules can change the configuration of their 
order parameter (e.g. from uniaxial to biaxial) and the EL equations cannot 
be used as a faithful model, which is rather given by the LdG theory based 
on the alignment tensor Q. Several dynamical fluid models for the evolution 
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of the alignment tensor Q are found in the hterature [1, 28, 29]. In this 
paper we shall show how the ideal Qiang-Sheng (QS) model [29] for the LdG 
tensor order parameter also allows for helicity conservation, in analogy to EL 
theory. 

The ideal QS model reads as 

dtU + u-Vu=-di(^-^-VQi,^ -Vp, V-n = (18) 
dtQ + u -VQ = J-^ P (19) 
dtP + u-VP = -% + - AI (20) 

where P is conjugate to Q and I is the identity matrix, while A is a La- 
grange multiplier arising from the condition Tr Q = 0. Here the free energy 
T{Q, VQ) contains the Landau-de Gennes free energy [9] as well as inter- 
action terms with external fields. Notice that the molecular field dJ^/dQ — 
di{dJ^/dQ^i) is always symmetric, so that P is also symmetric at all times. 
The circulation vector Cqs for the above system is defined by 

Cqs := It + Pjj VQjj . 

Indeed, a direct verification shows that the above vector satisfies equation 
(8), that is (cf appendix A. 2) 

dtC^^s + V (w ■ Cqs) - n X V X Cqs = -V0 , (21) 

with 

Thus, the Euler-like equation 

dtUcis + V X (u X cJqs) = (23) 
holds for the modified vorticity oJqs = V x Cqs. The circulation theorem 

^^^^(ii + P,,VQ,,) ■dx = 0, (24) 
and the helicity conservation (for lJqs and u both tangent to the boundary) 



d_ 

di 



Cqs ■ cUqs d^x = 
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are a natural consequence of the Euler-like equation (23) for the QS model of 
LdG theory. Again Ertel's commutation relation (13) for uJqs follows easily 
from equation (23). Moreover, vortex structures similar to those appearing 
in EL theory also exist in the LdG formulation. 

At this point, one can ask about other LdG formulations and in partic- 
ular one wonders whether the latter also exhibit vortex structures and con- 
servation of hydrodynamic helicity. Among the LdG formulations of liquid 
crystal dynamics, the one by Beris and Edwards [1] is probably among the 
most common, although it is not known to possess helicity conservation. In 
particular, this theory treats the order parameter as a "conformation tensor 
field" , so that the symmetric matrix Q is replaced by a symmetric tensor field 
on physical space. This deep geometric difference is probably responsible for 
the absence of the hydrodynamic helicity in the Beris-Edwards formulation. 
In this sense, the peculiarity of the QS model for the LdG tensor dynamics 
resides in exhibiting helicity conservation and its associated vorticity dynam- 
ics. These quantities both involve coupling between velocity and structure 
fields. The next section shows how this is actually a situation common to all 
fluid systems exhibiting rotational symmetry breaking. 



4 Completely broken symmetries 

The tensor order parameter Q arises as usual from the broken rotational 
symmetry that is typical of liquid crystal materials. When this rotational 
symmetry is fully broken, one needs to account for the dynamics of the whole 
particle orientation, which is determined by an orthogonal matrix O (such 
that = O^). This is a situation occurring, for example, in spin glass 
dynamics [12, 5, 22]. Eringen's wryness tensor (here denoted by n) is written 
in terms of O as [11] 

4 = \ ^"""^ Omk d^Onk = ^ e™"' diOnk Ofc„\ , (sum over repeated indexes) 

(25) 

(In this section we suppress labels such as EL or QS, in order to better adapt 
to the tensor index notation.) Although orthogonal matrices are difficult to 
work with analytically and the use of quaternions could be preferable, the 
correspondence between quaternions and rotation matrices is not unique. 
Thus, following Eringen's work [11], we identify molecule orientations with 
orthogonal matrices. 
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It is the purpose of this section to show how equation (10) is not pecuhar 
of nematic hquid crystals. Rather, equations of this form are pecuhar of all 
systems with broken rotational symmetries. In particular, equation (10) also 
holds in the case of complete symmetry breaking, for a vorticity variable UJ 
depending on the fluid velocity w, on the full particle orientation O and on 
the angular momentum vector 

^rmri Omk ^ nk ^rmn ^ nk ^ km (^^) 

where ^ is the variable conjugate to O. In the well-known spin glass theory of 
Halperin and Saslow [1-3], the rotation matrix is small and thus it is replaced 
by its infinitesimal rotation angle 6, where exp(0) = O. Then, 6 and cr 
become canonically conjugate variables, as shown in [12]. Here we consider 
the whole matrix O to account for arbitrary rotations. 

In the case of fully broken symmetry, the (incompressible) equations of 
motion can be written for an arbitrary energy density £{cr^ O) as [22] 

d£ d£ 

dtu + u-Vu = -ffrV— + 7^7:— V0„„ - Vp (27) 

dar dOmn 

dtOrnn + ' VO^n = —^mkj-Q^Okn (29) 

Notice that the validity of the above set of equations is completely general. 
Indeed, the above system is derived in [22] in a general fashion, under the 
only hypothesis that the broken symmetry group is 5*0(3). More general 
broken symmetries can be certainly treated in the same way, although this 
paper focuses only on rotational symmetries. 

In this context, equation (7) generalizes immediately by considering the 
wryness tensor in (25). Thus, the new circulation vector is defined by 

C := U + (J ■ K = U + ]^ Ernns Omk'^Onk , (30) 

where u, cr and O satisfy equations (27), (28) and (29) . 

At this point it is natural to ask whether the new C satisfies equation 
(8). Remarkably, a positive answer again arises from a direct calculation by 
using the ordinary properties of the Levi-Civita symbol. One obtains (see 
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appendix A. 3) 



dtC + Viu-C)-uxV xC = -V {^p-^ \ufj . (31) 

Consequently, the Euler-like equation (10) holds also in this case when the 
rotational symmetry is completely broken. Explicitly one writes 

(9t(V X C) + V X (w X V X C) = 0. (32) 

In turn, equations (32) and (30) imply the circulation law 

^ ^ + ^ ^mns OS Omk^Onk^ ■ dx = , (33) 

and the helicity conservation 

d_ 
It 



C ■ w d^x = 

V 



with ZU = V X C Thus, the existence of vortex configurations is independent 
of the type of symmetry breaking characterizing the fiuid. Therefore, such 
vortices may exist in liquid crystals independently of the choice of order 
parameter. However, one should also emphasize that the energy conserving 
assumption may fail in several situations and one would then be forced to 
consider viscosity effects. Moreover, polymeric liquid crystals do not seem to 
fit easily into the present framework; rather their description requires other 
liquid crystal theories such as the celebrated Doi theory [10]. 



5 Geometric origin of the helicity invariant 

In the previous sections, the helicity and vorticity of various systems with 
broken symmetry have been presented. However, the explicit formulation of 
these results still lack some more justification that can be found in the deep 
geometric nature of these systems, as it was emphasized in [7]. This section 
aims to give a brief overview of the geometric setting of the liquid crystal 
equations that eventually leads to the explicit formulation of their helicity 
and vorticity. This will show how these quantities can be found without any 
of the calculation presented in the Appendix, by simply relying on geometric 
symmetry concepts. The reader is also addressed to [26, 14]. 
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As our starting point, we write the total Poisson bracket for a general 
(incompressible) fluid system with broken symmetry, involving an order pa- 
rameter space M. In this case, the dynamical variables consist of the fluid 
momentum m(x), the order parameter state Q(x) G M and its conjugate 
variable 'P(x), so that (Q(x),'P(x)) G T*M. For simplicity, we restrict to 
consider the case when M is a matrix vector space. The total Poisson bracket 
reads as 

where the angle bracket denotes the pairing 

The above bracket is derived from the relabeling symmetry that character- 
izes all fluid systems. In particular, this bracket characterizes all Hamil- 
tonian fluid systems with broken symmetry. The relabeling symmetry car- 
ried by the fluid emerges mathematically as an invariance property of the 
Hamiltonian functional ^ : T*Diff(M3) x T*C^{R^, M) under the dif- 
feomorphism group Diff(M^) of smooth invertible maps. Here the notation 
C°°(M^,M) stands for the space of M-valued scalar functions, i.e. the space 
of order parameter flelds. The reduction process induces a reduced Hamil- 
tonian H : X*(M3) X T*C~(M^ M) M, where X*(R^) denotes the space of 
differential one- forms, i.e. the space of fluid momentum vectors m(x). This 
process leading to the reduced Hamiltonian H = H{m, Q,V) is widely ex- 
plained in [26, 14, 23]. 

Each term in the above Poisson bracket possesses a precise geometric 
meaning. While the flrst term coincides with the Poisson bracket for ordi- 
nary fluids, the second term is the canonical bracket for the order parameter 
held Q(x) and its conjugated momentum 'P(x). Moreover, the whole second 
line contains the two terms arising from the action of the relabeling symme- 
try group Diff(M^) on the canonical order parameter variables (Q(x),P(x)). 
Poisson brackets of this form were applied in different contexts, from elec- 
tromagnetic charged fluids [16, 15] to superfluid dynamics [18], and even to 
superfluid plasmas [19]. 
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At this point, upon following the Hamiltonian version of Noether's theo- 
rem (see [26]), one can construct the total momentum 



C = m + J{Q,V) 



where J{Q,V) is the (cotangent- lift) momentum map of components 



The geometric meaning of this momentum shift by a momentum map is best 
explained in [23]. In Lie derivative notation, the dynamics of the circulation 
quantity reads as 



which then arises naturally as the circulation conservation determined by the 
relabeling symmetry of the system. The explicit proof of circulation theorems 
of this kind can be found in many works in geometric fluid dynamics; see 
[20] for a modern reference. After recalling that for incompressible flows 
m = it is easy to recognize that replacing M by the space SymQ(3) of 
traceless symmetric matrices transforms the relation (36) exactly into the 
relation (21), which then produces the results in Section 3. Moreover, the 
corresponding vorticity relation for cJ = dC is easily obtained by taking 
the exterior differential of equation (36) and recalling that this operation 
commutes with Lie derivative. Then, one obtains {dt + £u) dC = 0. The 
form of the helicity is also easily derived from the above arguments, upon 
recalling an old result in [23]. In particular, if 'H{m) denotes ordinary Euler's 
helicity, then 'H(C) is a Casimir invariant of the Poisson bracket (34). Notice 
that all the above relations hold for an arbitrary manifold M other than 
a matrix space. This only requires using the appropriate pairing between 
vectors and co-vectors. 

So far, we only used the cotangent-lift momentum map, which can be 
found for all the cases when the dynamics involve conjugate variables in a 
cotangent bundle T*M. However, this does not appear to be the case for 




(36) 



thereby yielding Noether's conservation relation 
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the discussion in Section 4, where M = SO (3) and Q = O. This apparent 
contradiction is easily solved by noticing that 

TT{V^diO) = TT{{VO-YdiOO~^) . 

Then, upon denoting a = V0~^ and = diOO~^, the usual isomorphism 
between antisymmetric matrices in the Lie algebra so (3) and vectors in 
yields the term cr ■ k in the circulation quantity (30). Then, upon repeating 
the same steps as above, the momentum map TiiV^diO) = cr ■ k returns 
exactly the same results as in Section 4. 

The case of nematic liquid crystals treated in Section 2 can be also ob- 
tained by a direct computation, upon setting M = 5*^/^2, which is the 
director space. Upon denoting tt = JDtn the corresponding conjugate vari- 
able, it is easy to see that Vn ■ tz = cr ■ 7el. However, the geometric meaning 
of this simple step requires more basis that can be found in [7], where this 
last relation is justified by Lagrangian reduction. 

At this point, it is clear that the above arguments ensure the results in this 
paper without any need for further discussion. Nevertheless, the Appendix 
gives explicit proofs that can be followed without previous knowledge in 
geometric mechanics. 

6 Conclusions 

This paper provided explicit expressions for the helicity conservation in liq- 
uid crystals, in both EL and LdG theories. This conservation arises from 
an Euler-like equation that allows for singular vortex structures in any di- 
mension. Some of the ideal fluid properties were extended to liquid crystal 
flows, e.g. Ertel's commutation relation. These results were also shown to 
hold for molecules of arbitrary shapes, by considering fully broken rotational 
symmetries occurring in some spin glass dynamics. All of the results were 
eventually justified by geometric symmetry arguments. 

The energy- Casimir method can then be applied to study nonlinear sta- 
bility properties of these systems, see [21] for several examples of how this 
method applies to many types of fluids. This can be used, for example, to 
explore the coupled macro- and micro-motion of the stationary (generalized 
Beltrami) solutions. While the 3D stability analysis is limited by the fact 
that the helicity is the only Casimir invariant, the 2D case is much reacher 
because the whole family (1) of Casimir invariants becomes available. 
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One more remark concerns physically observable effects. More particu- 
larly, one wonders how conservation of total circulation causes observable 
effects. Even more, one would like to observe these effects in a particular 
experiment. A simple technique that could be used to this purpose is the use 
of external electric fields that drive the order-parameter variables, thereby 
generating fluid circulation by conservation of the total circulation. Then, if 
one applies an external field to a trivial motionless liquid crystal, the director 
alignment caused by the field would result in the generation of fluid motion. 

Other physical questions also arise about the nature of vortex solutions, 
which evidently represent much more that simply disclinations dragged around 
by a smooth flow. It is possible that these solutions share many analogies 
with superfiuid vortices in He^-A, whose order parameter space is again the 
whole group 5*0 (3). However, the nature of these singular solutions is left 
open for future investigations, together with their stability properties. 

Acknowledgements. We are indebted with Darryl Holm for his keen re- 
marks about the relation between the modified vorticity and helicity conser- 
vation in Ericksen-Leslie theory. Some of this work was carried out while 
visiting him at Imperial College London. 



A Appendix 

A.l Derivation of equations (8) and (10) 

Upon using the notation £u for the Lie derivative with respect to the velocity 
vector field u [26], we can rewrite equations (4)- (5) as 

^ + i^u^ w = Vn ■ h - V (^p + F - ^ |u|'^ (37) 

^ + i:^ j o- = h X n , + £^ j n = J" V x n. (38) 
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Then, one simply calculates 



(^^ + £u^ C^L = + (tt + cr^L ■ n X Vn) 



= Vn - h — V^J9 + F — + h X n ■ n X Vn 

- J~Vel • Vn X (ctel X n) + J^Vel ■ n x (Vctel x n) 
+ J^Vel ■ n X (cTel X Vn) . 

At this point, standard vector identities yield 

h X n ■ n X Vn = (h ■ n) (Vn ■ n) — (Vn ■ h) (n ■ n) 
= -Vn ■ h 

Vn X (ctel X n) = (Vn ■ n) ctel - (Vn ■ ctel) n 

= - (Vn ■ CTel) n 
n X (Vctel X n) = (n ■ n) Vctel - (Vctel ■ n) n 

= VcTel - (VcTel ■ n) n 
n X {(Tel X Vn) = (Vn ■ n) ctel - (ctel ■ n) Vn 

= 0, 

where we have made use of the relations |np = 1 and ctel ■ n = 0. Therefore 
equation (8) follows directly from 

The equation (10) follows by taking the curl of the above equation, upon re- 
calling that the curl is given by the exterior differential, so that d(CEL ■ dx) = 
(V X Cel) ■ dS. Since the differential always commutes with the Lie derivative 
[26], equation (10) follows immediately. It is also easy to see that equation 
(15) arises by calculating {dt + £u) (<^el ■ ^el) = ''^el ■ V0. 
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A. 2 Derivation of equations (21) and (23) 

Upon introducing the Lie derivative notation, the equations (18)-(19)-(20) 
read as 

(^^ + £u^u = hijVQij - V (^p + - ^\u\^^ , \/-u = (39) 
^^ + £u^Q = J^' P (40) 
^ + £^]p = -h-\l (41) 



9t 

where we have denoted the molecular field by 

dQ ""'dQ, 

Then, one simply calculates 



^ + £)c -(^ + £ 



which becomes 



^ + £^^ = -v(^ + F + \5,,Q,, - ]^\u\^ - ^P.,P., 
Finally, taking the curl of the above equation returns (23). 

A. 3 Derivation of equation (31) 

Upon using the Lie derivative notation, equations (27)- (28)- (29) may be writ- 
ten as 

d \ d£ d£ / 1 \ 

di^ = -^^^^ + - ^ - 2 J ^^2) 

^ + ^.J-. = %.^-.^-O..^J (43) 

' d \ dS 

TTT ~^ £u 1 Ojjiji ^ink'jT{ ^fcn (44) 

dt J 'da, 
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d \ f d 



SO that, upon denoting (j) = p — \u\'^/2, one computes 

d£ d£ 
= - cTrV- h — V0„„ - V0 

1 / n -^\n 

~l~ „^mns I ^sji^i r% ^sji^ih I ^mk^^nk 

1 9^ 
2 

1 „ fv7^^ V7^ 

2 \ 9(Tj ao-j 
At this point we observe that, since OmkOhk = OmkO^^l = Smh, then 

(7(Tr 2 (7(Tj acr^ 2 oaj 

= -CTr Vt^ h Osj OsV^r— = 0. 

Moreover, the sum of all terms in dS/dcr can be written as 
d£ 

C"s i^^ijs ^mni Omk^ Ofik ^ims ^inj Onk^Omk £mis ^inj^ Onk Omk) 

88 

CTs T; ^ijs ^imn Omk^ Onk 

dS 

(^s i^^jm ^sn ^jn ^sm) Omk^ Onk 

d£ 

= -^Sg^^{VO.kO^}+OskVO,}) 

d£ 

= -^^^v(o.,o,7)=o. 
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In addition, we calculate 



d£ 1 d£ 

dOmn 2 dOjh 

dS 1 dS 

1 

— V0„„ + - (Oi/, 0,fe VO.fc - Oik VO.fc) — — 
90™„ 2 dOjh 

8F 1 (9/^ 



"90„„ ' 2 ^ '''' '^^^ ■'^^ dOjh 

d£ 1 (9£ 

OS „^ „^ 

V0™„ - VOjh^^ = 0. 



dOmn dOjh 

Thus, we have proved the relation 

d \ 

- + £^]C = -V<j), 

whose curl yields the corresponding Euler-like equation (32), thereby recover- 
ing the corresponding helicity conservation. Notice that the above result also 
holds in the case of explicit dependence of the free energy S on the gradient 
VO of the orientational order parameter. 
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